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Abstract 

We show the possibility of quasiperiodic waves at the onset of 
thermal convection in a thin horizontal layer of slowly rotating zero- 
Prandtl number Boussinesq fluid confined between stress-free conduct- 
ing boundaries. Two independent frequencies emerge due to an inter- 
action between a stationary instability and a self-tuned wavy insta- 
bility in presence of coriolis force, if Taylor number is raised above a 
critical value. Constructing a dynamical system for the hydrodynam- 
ical problem, the competition between the interacting instabilities is 
analyzed. The forward bifurcation from the conductive state is self- 
tuned. 

PACS number(s):47.20.Ky, 47.27.Te 



Thermal convection in Boussinesq fluids in the limit of vanishing Prandtl 
number P(= v/k) ([p|-||12||) is of interest for astrophysical problems (P ~ 
10~ 8 ) as well as for liquid metals (P ~ 10~ 2 — 10~ 3 ). The theoretical study, in 
particular with stress-free boundary conditions, in the limit of large thermal 
diffusitivity k has been considered subtle for a long time because the linearly 
unstable two-dimensional (2D) rolls become exact nonlinear solution. The 
nonlinearity v ■ W due to the advection of temperature fluctuation 9 by the 
velocity field v might be negligible IJ. The nonlinearity v • Vv, due to the 
self interaction of velocity field, does not contribute to saturation for straight 
(2D) rolls. This led to the speculation that the zero P limit might involve 
a singular limit problem similar to the one with infinite Reynold number in 
incompressible fluid dynamics. However, the recent three-dimensional (3D) 
direct numerical simulation (DNS) of zero-Prandtl number Boussinesq equa- 
tions by Thual [II] showed the saturation instead of indefinite growth of the 
solution even with stress-free boundary conditions. He also compared the re- 
sults of zero P equations with that of the full Oberbeck-Boussinesq equations 
in the asymptotic limit of vanishing P, and found complete agreement in two 
cases. The saturation of growing 2D rolls at the onset of convection occurs 
by generation of self-tuned 3D waves, the mechanism of which was explained 
in a simple model fl2 |. The results of this model, in its validity range, agreed 



well with that of DNS just above the convective instability. The new nonlocal 
instability at the onset occurs purely due to nonlinear effects, while linear 
equations predict the stationary instability [T3| . 

We present, in this article, a dynamical system constructed for thermal 
convection in zero Prandtl number Boussinesq fluid, confined between stress- 
free conducting flat boundaries, and subjected to a slow rotation about the 
vertical axis. We then investigate numerically the system to study the effect 
of coriolis force on the onset of convection. We show that the convection 
sets in as quasiperiodic waves at the onset of convective instability for Taylor 
number T above a critical value T c , although the principle of exchange of 
stability is valid according to the linear theory |13| for these Taylor numbers. 
The generation of two independent frequencies is the result of an interac- 
tion between a stationary instability and a self-tuned wavy perturbations in 
presence of coriolis force. This an example of a new self-tuned forward bi- 
furcation. For the values of Taylor number below T c , the convection sets in 
as self-tuned wavy instability as is the case in absence of rotation. However, 
the model shows a possibility transition from one wavy instability to another 
through a narrow window of period-doubling instability. 
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We consider a thin layer of a Boussinesq fluid of infinite horizontal exten- 
sion subjected to a uniform adverse temperature gradient (3 across the fluid 
layer, and a rigid body rotation with an angular velocity Q about the vertical 
axis. The fluid is assumed to have uniform values of the kinematic viscosity 
v and the thermal difusitivity k. The basic state is the conductive state with 
no fluid motion in the rotating frame of reference. The convective flow, in 
the limit of zero-Prandtl number, is then described by the following system 
of dimesionless hydrodynamic equations, 

d t (V 2 v 3 ) = \7 4 v 3 + RV 2 H 9 - VTd z u 3 

- e 3 -V x [(a; • V)v - (v • V)w] , (1) 

d t u 3 = V 2 u 3 + Vfd z v 3 

+ [(u ■ V)v 3 - (v • V)w 3 ] , (2) 

V 2 £ = -v 3 , (3) 

where v(x,y,z,t) = (vi,v 2 ,v 3 ) is the velocity field, 8(x,y,z,t) the deviation 
in temperature field from steady conduction profile, u = (ui, cu 2 , u 3 ) = V x v 
the vorticity field of the fluid. In the above, length scales are made dimen- 
sionless by the thickness d of the fluid layer, time by the viscous time scale ^, 
and the temperature field by ((3(£)~- Rayleigh number R = a9 ^ d and Tay- 
lor number T = 4d4 ^ 2 are the two dimensionless external parameters. The 
unit vector e 3 is directed vertically upward. We impose periodic boundary 
conditions in horizontal plane. This introduces two fundamental wave num- 
bers k along x-axis and q along y-axis. The stress-free boundary conditions 
imply d z v\ = d z V2 = v 3 — at z — 0, 1. Thermally conducting horizontal 
boundaries mean 9 — at z = 0,1. 

The hydrodynamical equations (|l|-|3D are the same as those derived by 
Chandrasekhar [13|. We have nondimesionlized them and considered the 



case of zero P. We have also eleminited the pressure term from Navier- 
Stokes equations by taking curl twice and using the incompressibility condi- 
tion (V ■ v = 0). The conclusions derived from the linearized version of the 
equations remain unchanged even in the present case. Following the argu- 



ments of Chandrasekhar |T3| , one arrives at the conclusion that the principle 



of exchange of stability is valid even in the limit discussed here. The critical 



value of Rayleigh number now reads R C (T) 



and crit- 



ical wave number k c (T) = tx^J {l\ + l 2 — \) now depend on Taylor number 
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T [O. In the above, 



and x = k^/ir 



i {i + S ± y(f + S) ~~ i 

In absence of rotation = fc c (T = 0) = n/y/2 and = -R C (T = 0) 
27tt 4 /4. 

The 2D rolls are not exact solutions of nonlinear hydrodynamic system 
with rotation as is the case in zero-Prandtl number convection in absence of 
rotation. Nevertheless, the growing 2D rolls are not saturated just above the 
onset of convective instability. The saturation occurs only because of nonlin- 
ear interaction of 2D rolls with 3D wavy perturbations. To understand the 
nonlinear behavior close to the convective instability, we construct a con- 
sistent minimal- mode model using Galerkin technique [R]. We expand the 
vertical velocity and the vertical vorticity 003 in Fourier series compatible 
with the stress-free boundary conditions and conducting thermal boundary 
conditions. As the DNS, in absence of rotation, showed standing patterns 
[ 12 1 instead of traveling patterns, we expect similar behavior at least for small 



rotation rates. Therefore, we expand the fields with real Fourier coefficients. 
This lead to the following expansion for the vertical velocity and the vertical 
vorticity in a minimum-mode model. 

v 3 (x,y,z,t) = W W i(t) cosk c xsmnz 

+ Wm (t) cos k c x cos qy sin nz 

+ Wni(t) smk c xsm qy sin ttz 

+ W012 (t) cos qy sin 2nz + . . . (4) 

u> 3 (x, y, z, t) = Cioi(^) cos k c x cos7rz 

+ C010 (*) cos qy + Clio (*) sin k c% sin qy 

+ Cm(^) cos A; c a; cos g|/ cos 7T2; 

+ (012 (i) cos qy cos 2nz + C2oo(^) cos 2k c x 

+ C210 (t) cos 2 k c x cos qy 

+ C210W sin 2k c x sinqy + ... (5) 

The mode selection is quite systematic. As rotation couples the vertical 
velocity and the vertical vorticity linearly, we have selected the mode Cioi- 
The mode £010 is essential to saturate zero Prandtl number convection via 
wavy instability. All other modes appear through the nonlinear interaction of 
these vorticity modes with the critical velocity mode Wioi- As the vorticity 
field is very crucial for saturation in the limit of vanishing Prandtl number, 
all relevant second harmonics are retained for vertical vorticity. All relevant 
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harmonics of the vertical velocity field, consistent with the selection of the 
vertical vorticity, are also retained. Other higher order modes may be re- 
quired as Rayleigh number is raised further. As we are interested to capture 
essential nonlinear interaction between competing instabilities just above the 
onset of convection, these modes are essential. The solenoidal character of 
the velocity and the vorticity fields yield horizontal components of the ve- 
locity and the vorticity fields. The thermal fluctuation 9 is captured from 
Eq. [3[ Projecting the hydrodynamic equations ([I] - ||) on above modes, we 
arrive at a twelve-dimensional dynamical system [|i~5|j . 

We now investigate the solutions of the dynamical system by performing 
numerical integration of the model using standard fourth order Runge-Kutta 
as well as Bulirsh-Stoer schemes. By choosing a value for T, we set k c (T). 
We then choose a value for q. We have tried with different values of the 
wavenumber q of the perturbations and got qualitatively similar results ex- 
cept when the ratio q/k c is close to unity. We present here all the results 
for the case q/k c (T) = 0.4. Initial values for all the twelve modes are chosen 
randomly, and integration is done for a fixed value of Rayleigh number R. 
We then repeat the process by increasing the value of R in small steps. We 
have also tried various initial conditions. The results of all the numerical 
integrations remain the same for the same values of all the relevant parame- 
ters. In absence of rotation (T = 0), only six modes are excited. This model 
then reproduces the results of the model |12| of zero P convection without 
rotation. In presence of rotation all twelve modes are excited as it should in 
a consistent model. 

Figure 1 gives the stability boundaries of various possible solutions, in 
the parameter space R — T, computed from the model dynamical system. 



The lowest line in this figure shows linear dependence |13[ of R c , the critical 
Rayleigh number, on Taylor number T for the onset of stationary convection 
in zero Prandtl number Boussinesq fluid. The onset of overstability for the 
case of vanishing P for the Taylor numbers considered here is much above, 
and is not shown in the figure. 

As Rayleigh number is raised above its critical value R C (T) for various 
values of Taylor number below T = 6.0, conduction state becomes unstable 
via stationary bifurcation. However, 2D rolls with broken mirror symme- 
try |16] does not saturate until wavy perturbations interact with them. This 
saturates the growing rolls at finite amplitude even at the onset. The wavy 
perturbations are automatically generated when the amplitude of the 2D roll 
mode becomes large enough. The self-tuned waves consume the energy of 
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2D rolls and stop the unbounded growth of latter. This is precisely what 
happens when there is no rotation and 2D rolls have mirror symmetry. 

The three rows of Fig. 2 show projections of the phase diagram, starting 
from left, in £ 10 i — Wioij Coio ~ Wioi> and Wm — planes for various 
Rayleigh numbers (increasing downward) for fixed value of T and q. 

As Rayleigh number is increased slowly, the solution changes from one 
wavy solution to another through a thin regime showing period doubling 
solutions (see the middle row of Fig. 2). The first wavy solution SW1 has 
2D mode W\oi with non-zero mean as in the absence of rotation in zero P 
convection ||12|| , while the second wavy solution SW2 has 2D mode with zero 
mean. As Rayleigh number is increased, the exchange of energy from 2D 
modes to waves increases. The larger amplitude variation of vorticity modes 
is at the cost of energy of 2D rolls. This is well known feature in the case of 
oscillatory instability. 

As T is raised further, the rotation facilitates easily the exchange of 
more energy from the 2D roll mode Wioio to the vertical vorticity mode 
C101 through linear coupling. We observe an interesting behavior for T > 6 
(see Fig. 1). The conduction state becomes unstable via stationary insta- 
bility |[3j but the final state just above onset is quasiperiodic waves |17| . 
Figure 3 shows the variation of various modes with time. The amplitudes of 
all the modes begin modulating at the same frequency. However, the ratio 
of the frequency of wavy motion and that of amplitude modulation not an 
integer. The Fourier transform of these modes shows two independent fre- 
quencies. The frequency of amplitude modulation is much smaller compared 
to that for wavy motion. The sharp decrease of the amplitudes of higher 
order modes confirms the fast convergence of the expansion. The model, 
therefore, represents accurately the scenario close to the instability onset. 
Figure 4 shows the projections of phase space trajectories in various planes. 
It clearly describes the quasiperiodicity of the convective flow. The trajec- 
tories are confined in twelve dimensional torous in the phase space. The 
quasiperiodic behavior originates due to the nonlinear interaction among the 
2D velocity mode W W i, the 2D vorticity mode ( W i excited by rotation, and 
the wavy vorticity mode £oio- Figure 5 reveals some interesting details of 
the time dependence of of convective patterns. The complex textures of the 
quasiperiodic patterns are shown over a period of wavy motion, which is 
much faster than the amplitude modulation. The halves of a period of wavy 
motion are quite asymmetric. The textures of the pattern at different times 
are never the same due to quasiperiodicity. 
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We have presented in this work a simple dynamical system, which de- 
scribes the phenomenon of thermal convection in rotating Boussinesq fluid 
of zero Prandtl number very close to the onset. For Taylor number above 
a critical value T c , quasiperiodic waves are observed at the instability onset. 
For very values of Taylor number below T c , the coriolis force causes one wavy 
instability to another through period doubling instability. We have shown 
that convection might be possible as quasiperiodic waves, even if the prin- 
ciple of exchange of stability is valid according to linearized hydrodynamical 
system. The saturation to quasi-periodic convectie state is self-tuned purely 
due to the nonlinear effects. This is an example of new self-tuned bifurcation 
scenario from the conduction state to unsteady convective state. The model 
presented would also be useful to study an interesting possibility of transition 
from a state of rest to quasi-periodic chaos (T^j at the primary instability. 

Acknowledgements: This work was sponsored by DST, Govt, of India 
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Figure Caption 



Fig. 1: Stability boundaries in the parameter space R — T just above the 
onset of convection. The conduction state is stable below the lowest straight 
line, which shows critical value R c of Rayleigh number as a function of Taylor 
number T (q/ K c = 0.4). The region marked QP shows quasiperiodic waves. 
The regions denoted by SW1 and SW2 show two regimes of wavy solutions 
separated by a thin region where period-doubling is observed. 

Fig. 2: Phase portrait for various values of Rayleigh number R (T = 2.0 
and q/k c = 0.4). The top row, starting from left, shows the plots of modes 
Cioi) Coio? and Win with respect to W W i for R = 663.0. The middle and the 
bottom rows show the same plots for R = 664.95 and R = 668 respectively. 
The top and the bottom rows two different wavy regimes SW1 and SW2 
respectively, while the middle row represents period doubling corresponding 
to the narrow regime between SW1 and SW2 in Figure 1. 

Fig. 3: Time variation of all the modes for T = 10.0, q/k c = 0.4, and 
R = 678.0 long after all transients are died out . The critical Rayleigh num- 
ber R C (T = 10) = 677.0768. Starting from left, the top row shows variation 
of Wioi, Wiii, Win, and Wqu with time. The middle row shows variation 
of Cioi, Coio, Cm, C200, and the bottom row shows variation of C012, (210, C210, 
and Clio- 
Fig. 4: Phase space portraits showing quasiperiodic motion for T = 10.0, 
q/k c = 0.4, and R = 677.08. Starting clockwise from the left top, they show 
the projections of of the phase space in W W i — W\\i, C101 — W^ioi? C210 — W W i, 
and Wqyi — Wiqi planes. Some modes like C210 and Wou have faster time 
period (i.e., period of wavy oscillation) double that of W101. 

Fig. 5: Contour plots (T = 7.0, q/k c = 0.4, R = 671.30) at z = 0.25. 
Various textures of the pattern is shown over a period of faster time scale 
to at the equal time interval of to/8. The sequence of time evolution of the 
pattern-texture is shown from left to right in each row starting from the top 
row. 
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